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A STUDY OF THE APPLICATION OF POWER-SPECTEAL METHODS OF GENERALIZED 
HARMONIC ANALYSIS TO GUST LOADS ON AIRPLANES ^ 

By Habht Press and Bernard Maiblskt 


SUMMARY 

The applicability of some resvlts from the theory of general- 
ized harmonic analysis (or power-spectral analysis) to the 
analysis of gust loads on airplanes in continuous rough air is 
examined. The general relations for linear systems between 
power spectrums of a random input disturbance and an output 
response are used to relate the spectrum of airplane load in 
rough air to the spectrum of atmospheric gust velocity. The 
power spectrum of loads is shown to provide a measure of the 
load intensity in terms of the standard deviation {root mean 
square) of the load distribution for an airplane in flight through 
continuous rough air. For the case of a load output having a 
normal distribution, which appears from experimental evidence 
to apply to homogeneous rough air, the standard deviation is 
shown to describe the probability distribution of loads or the 
proportion of total time that the load h&s given values. Thus, 
for an airplane in flight through homogeneous rough air, the 
probability distribution of loads may be determined from a 
power-spectral analysis. 

In order to illustrate the application of power-spectral 
analysis to gust-load analysis a/nd to obtain an insight into the 
relations between loads and airplane gust-response character- 
istics, two selected series of calculations are presented. In the 
first series, the standard deviations of loads in continuous rough 
air described by an assumed power spectrum are calculated for 
systematic variations in the frequency and damping character- 
istics of the airplane response to a step-gust input. The results 
obtained indicate that the loads in rough air are particularly 
sensitive to variations in the damping characteristics of the 
oscillatory response to a step gust and largely independent of 
variations in the frequency. In the second application, the 
standard deviation of loads is calculated for selected variations 
of each of several airplane geometric and aerodynamic param- 
eters of an idealized and stable transport-type airplane. The 
standard deviations obtained are compared with results derived 
by conventional techniques of using the calculated peak response 
to an idealized and representative discrete gust. The results 
indicate that for stable conflgurations both methods of analysis 
yield results that are consistent to a first approximation. 


INTRODUCTION 

The study of gust loads on airplanes is a twofold problem 
requiring the adequate representation of the characteristics 
of atmospheric turbulence and the determination of the air- 
plane response (loads or motions) in rough air. These 
problems have been recognized since the inception, of gust- 
load research but because of the difficulties involved only 
li m ited approaches to the problems appeared practical. 
The methods that have been used are described and discussed 
in detail in reference 1. In general, the approach has been 
to use simplified airplane-response theory for the determina- 
tion of the characteristics of discrete gusts from airplane 
measurements of load. The “gusts” derived on thig basis 
are then used to calcffiate loads on other airplanes. Although 
these procedures appear reasonable for transferring loads to 
similar airplanes, as indicated in the reference, they are of 
questionable value for airplanes of widely different character- 
istics (such as, configurations and stability characteristics). 
These limitations have, however, not proved serious in the 
past since the transport airplanes which were primarily 
affected by gust standards appeared, in general, to follow 
conventional design. . Available data indicate that new 
transport airplanes experienced gust loads which were, in 
general, compatible with those predicted from past work. 

Trends in aeronautics toward higher speed and the develop- 
ment of missiles have served to introduce a widening range 
of unusual configurations and aircraft stability character- 
istics. Furthermore, the gust-load design requirements, 
which formerly were of concern for transport and bombei- 
airplanes only, appear to have become important for other 
aircraft as well. As a consequence of these developments, 
the need for more generally apphcahle techniques, both for 
the measurements of the characteristics of atmospheric gusts 
and for the calculation of the gust loads on new airplanes, 
has become more uigent. 

Developments in the theory of generalized harmonic 
analysis (ref. 2) appear adaptable for extending the analysis 
of gust loads beyond the discrete-gust case to the case of con- 
tinuous turbulence. Techniques from generalized harmonic 


> Supersedes NACA TN 2863, "A Study of the Application of Power-Spectral Methods of Generallicd Harmonic Analysis to Oust Loads on Airplanes” by Harry Press and Bernard 
Moiolsfey, 1063. 


371 



372 


EEPOBT 1172 — NATIONAL ADVISOBY COMMITTEE FOB AEBONATJTICS 


analysis involving the concept of power spectral density have 
been used for many years in diverse fields, such as in the 
study of random-noise problems in communications and iu 
the study of small-scale turbulence of wind tmmels^ The 
concept of the power of a random disturbance, which is 
fundamental to the present study, is defined by analogy to 
electrical power to be the time average of the square of the 
disturbance. The portion of the power arising from compo- 
nents having harmonic frequencies between w and is 

denoted as the power spectral density. For linear systems, 
the power-spectral-density fimctions of a random input dis- 
turbance and an output response are related through the 
frequency-response characteristics of the system. 

Attractive features of spectral analysis for the study of 
gust loads are the possibilities that: 

(1) Continuous turbulence can be described in analytic 
form by a power spectrmn rather than by discrete gusts. 

(2) The load response of airplanes to continuous rough air 
can be evaluated. 

(3) The desirable response characteristics of an airplane 
for minifniring gust effects in continuous rough air will be- 
come amendable to analysis. 

In view of the attractive features of power-spectral- 
density methods of analysis, an investigation of the applica- 
bility and implications of these techniques to gust-load 
analysis was undertaken and the results obtained are re- 
ported herein. In this reporj;, the concepts and relations of 
generalized harmonic analysis are defined and their method 
of application to the gust-load problem is indicated. The 
applicability of the normal probability distribution for the 
representation of the probability distribution of loads in 
continuous rough air is considered. For the case of a normal 
distribution of loads, the standard deviation (root mean 
square) of the load history defines the probabdity distribu- 
tion of loads. It is indicated that the standard deviation 
may be determined from load power spectrum; thus the re- 
sults of a power-spectral analysis permit the determination 
of the probability distribution of loads. Finally, in order 
to illustrate the apphcation of power-spectral analysis to 
gust-load calculations and to obtain an insight into the rela- 
tion between loads in continuous rough air and airplane 
characteristics, two applications are presented. In both 
applications, the power spectrum of atmospheric turbulence 
obtained from fiight measurements (ref. 2) is used to repre- 
sent the tmbulence input. The first application is intended 
to represent the effects on gust loads of variations in airplane 
dynamie longitudinal stability. The second application is 
intended to illustrate the effects on gust loads of variations 
of some geometric and aerodynamic parameters of an ideal- 
ized transport airplane. The indicated variations in load 
intensity are compared with those derived by conventional 
techniques of using the peak-load response to an idealized 
discrete gust. 


SYMBOLS 

A( ) response to unit step disturbance 

a, slope of tail lift curve per radian 
Oto slope of wing lift curve per radian 

b, ko parameters of equation (31) 

c airplane wing chord, ft 

Ct airplane tail chord, ft 

I pitching moment of inertia about center of gravity, 
slug-ft* 

V reduced frequency, oic/V, radians/chord 

Za, horizontal distance from center of gravity of airplane 
to wing aerodynamic center, ft 
It horizontal distance from center of gi-avity of airplane 
to tail aerodjmamic center, ft 
N nmnber of observations 

An acceleration increment, g 

P probabfiity 

P(z) cumulative probability distribution of z, equation (28) 
8 distance, chords 

S wing area, sq ft 

St tail area, sq ft 

t time, sec 

T arbitrary value of t, sec 

T( ) frequency-response function with respect to argu- 
ment, 6), Q, or V 

U vertical gust velocity, ft/sec 

V airplane true airspeed, ft/sec 

W airplane weight, lb 

W{ ) response to unit impulse disturbance 
X distance, ft ♦ 

X arbitrary value of x, ft 

y arbitrary random variable 

/(y) probability density function of variable y, equa- 
tion (23) 

y{ ) arbitrary function of arguments t, x, and s 
y^{t) average power of y{t), equation (4) 

z standardized variable. An/a 

az coefficient of skewness, equation (30b) 

04 coefficient of kurtosis, equation (30c) 

a standard deviation of specified random variable, 

equation (24b) 

p mass air density, slugs/cu ft 

T time displacements, sec 

B(t) autocorrelation function, equation (9) 

$( ) power-spectral-density function of an arbitrary dis- 
turbance with respect to w or 

#,( ) power-spectral-density function of a designated input 
4>o( ) power-spectral-densityfunctionofa designated output 

f(z) normal distribution with mean of 0 and standard 
deviation equal to 1, equation (27) 
a cireular frequency, radians/sec 

Q reduced frequency, osjV, i-adians/ft 

dt/da downwash factor 
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Subscripts: 

vfiax m aximum load response to a discrete gust 
basic basic airplane configuration 

A bar over a symbol designates the average value of the 
quantity. 

In this analysis, the use of several independent variables t, 
X, and 8 for an arbitrary disturbance y( ) and their associated 
frequency argiunents w, Q, and v has been found necessary. 
In order to designate that the several functions y( ), $( ), 
and T( ) depend upon their ai^uments a circumflex a and 
a tilde ~ have been used over the appropriate sets of functions 
in accordance with the scheme shown in the foUowmg table: 


Function 

Vartable 

t, sec 

X, ft 

*, chords 


»C0 

T(f<S) 

TV® 

A® 

V(x) 

Q 

J(Q) 

T(m) 

Aw 

A (I) 

»(») 

P 

i(p) 

f(M 



FroqucnoyHcsponse function 




ANALYSIS 


BASIC CONCEPTS AND HELATIONS OF POWER.SPECTHAL ANALYSIS 

The theory of generalized harmonic analysis is an out- 
gi-owth and generalization of harmonic analysis and is 
largely the work of Norbert Wiener (ref. 3). Accounts of 
the theory and method of appHcation are also given in 
references 2, 4, and 5. In order to orient the reader, a brief 
accoimt of the background, the basic concepts, and the 
relations is presented. 

The theory of harmonic analysis indicates that an arbi- 
trary periodic fimction can be represented by a Fourier 
series in the following manner: 

F(0=^+S (-4* cos noit+B^ sin nut) (1) 

n**l 

where 

2 

A„’=y I F(t) cos 7uat dt 
Fn=^J F (0 sin nut dt 

and T is equal to 2 t/ci) and is the period of the function. In 
order to apply this technique to nonperiodic phenomena, 
the hmit of equation (1) as T goes to infinity must be con- 
sidered. For this case of a nonperiodic function, the Fourier 
series takes the form of the Fourier integral and is given by 

[J“F(r)e-‘-dr]dco (2) 

If the second integinl on the right is denoted as 0{u) then 


equation (2) can be written m the reciprocal form 

F{{)=^^y{u)e*’^^du (3) 

where 

G{u)= r F{T)e-*“^dr 

%J —03 

The Fourier transform pair thus provides reciprocal relations 
between the time function F{t) and its frequency represen- 
tation 0{u). The quantity' Q{u)du gives the contiibution 
of those harmonic components of F{t) whose frequencies 
lie between u and u-\-du. 

A necessary condition for the apphcation of equation (3) 
is that the integrals involved be convergent. This condition 
acts as a severe limitation on the apphcabflity of the Fourier 
integral relations. In many problems, such as noise in an 
electric circuit or turbulence encoimtered by an airplane 
in flight, the disturbance is nonperiodic, persists for a 
relatively long period of time, and shows no tendency of 
dying out. In these cases the required integrals do not 
converge and, as a consequence, the frequency representation 
in terms of equations (2) and (3) is not possible directly. 

In order to develop a frequency representation which 
would be applicable to continuing disturbances, the theory 
of random processes makes use of the concept of a stationary 
random process. The characteristics of a stationary random 
process are described in detail in reference 4. Essentially, 
the assumption is that the underlying mechanism which 
gives rise to the disturbance does not change in time and 
that a statistical equilibrium exists. Thus, the statistical 
characteristics of the distribution ore invariant with time 
and statistical prediction becomes possible. For the case 
of a stationary random function of time y(jt), the mean 
square y‘{t) is defined by 

W^=^kr{y(f)fdt (4) 

T*— »CD ^ 0 

The mean square will usually exist and represent a measure 
of disturbance intensity. Since is a quadratic function 
of y{f), it has been termed the “average power” of y{f) in 
analogy to electrical power which is proportional to the square 
of the current. The function y{t) is considered to be com- 
posed of on infinite number of sinusoidal components with 
circular frequencies u, between 0 and oo. The portion of 
y^(f) arising from components having frequencies between 
u and w-f-dcj is denoted herein by 4>(w)dw. The function 
$(a>) has been called the power-spectral-density function in 
the hterature. From this definition, #(w) has the property 
that 

$(w)dco (5) 

The power-spectral-density function of a random variable 
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y(t) is generally defined in the following manner (see, for 
example, ref. 4) 


power-spectral-density fimction is consistent with the 
preceding definition of equation (6).- 


1 PT 

4>(a>)= lim — 7 K y{t)e~*"‘dt 

T-^a, ■KJ. Jo 


( 6 ) 


where the nota,tion || indicates the modulus of the complex 
quantity. If $(co) is defined in this manner, it has the prop- 
erty that its integral over the limits 0 to oo is equal to the 
power y^(f). That this expression for $(w) is consistent 
with the preceding discussion is seen to he plausible from the 
fact that, for each frequency, #(&j) is proportional to the 
square of the amplitude of the component of y(tf at that 
frequency. Thus, $(w)<Zco is a measure of the contribution 
of that frequency to y^(t). It should be mentioned that the 
definition equation (6) differs in minor detail from that used 
in reference 2 but agrees substantially with those used in 
references 4 and 6. 

A significant and useful relation for linear systems exists 
between the power-spectral-density function of a 

random input disturbance and the power-spectral-density 
function #o(w) of an associated output through the system 
frequency-response fimction T(u3). The system frequency- 
response fimction (or admittance) T(i<a) is defined such that 
T(io))e*“‘ is the system response to the sinusoidal input 
e‘"‘. In these terms, the relation between the power- 
spectral-density functions is given by 

^o(«) =$,(«) ir(t«)l* (7) 


For a given linear system, the function T(ico) may he con- 
veniently obtained from either the unit impulse . response 
or the unit step response, respectively, by means of the 
following relations: 


T(icJ)= J" W (t)e' 
00 

T('i<o)=iw J Al(Q< 


-<«< 



where TF(t) is the response to a unit impulse and A(0 is the 
response to a unit step. 

Equation (6) may be used to evaluate the power-spectral- 
density fimction from observed data. However, in practice, 
the power-spectral-density fimction may he determined more 
conveniently and less tediously by using a related function, 
the autocorrelation function R{j), defined by 


B(t) = lim ^ r y (t) y(t-f r) dt 

T-»eo X Jo 


(9) 


The autocorrelation function has the symmetrical properly 
and is reciprocally related to the power- 
spectral-densily function by the Fourier cosine transforma- 
tion in the following manner: 



$(o)) cos cor dco 


y 


2 T" 

$(co)=— I B(t) cos cot dr 
vjo 


( 10 ) 


Eeference 4, for example, shows that this definition of the 


some forms of POWER.SPECTEAL RELATIONS FOE GUST APPLICATIONS 

In this section, the method of application to gust loads on 
airplanes of the aforementioned concepts and relations is 
considered. Some available information is presented on the 
frequency-reqionse function T('ico) for the gust-load condi- 
tion and on the power-spectral-density functions for atmos- 
pheric turbulence. Finally, for convenience in ' gust-load 
applications of power-spectral methods, the functions and 
relations axe presented in terms of distances rather than time. 
Use is made of the distances x in feet . and a in chords and 
their associated frequency arguments £] and v. These 
changes of variables are subsequently shown to be particu- 
larly appropriate in the gust case because in these terms the 
power spectrums are independent of airplane fonvard speed. 

If atmospheric turbulence can be considered a stationary 
random process, then the basic requirement for the applica- 
tion of the foregoing power-spectral concepts and relations 
is satisfied. Evidence that this assumption is plausible 
under some conditions exists and is discussed subsequently. 
On the basis of this assumption, the turbulent-vertical- 
velocity distribution along a line in space can at an instant 
in time be considered to represent a stationary random func- 
tion of space consisting of an infinite number of harmonics 
of various frequencies or wave lengths. For this condition of 
spacial wave lengths, a natural unit for the associated fre- 
quencies would appear to be radians per unit distance or 
radians per foot. The description of the power-spectral- 
density fimction of atmospheric turbulence must thus be given 
basically in terms of units such as radians per foot. However, 
for an airplane in flight through rough air, consideration of 
the airplane as penetrating the gusts and experiencing the 
associated loads in terms of time is frequently convenient. 
Thus, the power-spectral-density functions of gust velocity 
and loads experienced by the airplane may be considered 
expressible in terms of the frequency aigument w in radians 
per second used in the preceding section. It is, therefore, 
permissible to express the relation between the power spec- 
trums of gust velocity and loads or normal acceleration in 
terms of equation (7) as 

i>o(co)=$<(w) |2’(ico)l* 


where $o(w) is the power-spectral-density function of airplane 
normal acceleration, $<( 0 ) is the power-spectral-density 
function of gust velocity experienced by the airplane, and 
T{w) is the airplane normal-acceleration response function 
for a sinusoidal gust velocity input. 

The present study is primarily concerned vdth vertical 
gust velocity inputs and acceleration increment Ati (load 
factor) outputs. The function $<(w) wiU have the dimensions 


of 


(ft/sec) ^ 


|T(’Uo)p will be given in 


\ft/sec/ 


Conse- 


radian/sec' 

(7^ 

quently, $»(&)) will have the dimensions of In 

this form, the power-spectral-density fimctions, having 


the dimensions of — — , j-efer to a particular airspeed, 
radian/sec • 

Before considering the representation of these functions in a 
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form mdependent of airplane speed, some remarks on the 
detennination of the frequency-response function for the 
gust-load condition and on available information concerning 
the power-spectral-density function for atmospheric turbu- 
ence appear appropriate. 

Frequency-response function. — ^Experimental methods for 
determining the frequency-response function of an auplane 
for a gust velocity disturbance have unfortunately not yet 
been developed. The frequency-response function for air- 
planes for a gust disturbance can, however, be estimated by 
theoretical methods by solution of the airplane equations of 
motion for a gust disturbance. Methods for the determina- 
tion of the frequency-response function for a linear system 
are described in chapter 2 of reference 4, for example. The 
calculation of the frequency-response function may some- 
times be more conveniently performed by first determining 
the airplane respopse to a step gust by methods such as those 
described in reference 6. The response of a step gust may 
then be used in equation (8) in order to determine the 
frequency-response fimction for a continuous sinusoidal gust 
input. 

Power spectrum of atmospheric turbulence. — The power- 
spectral-density functions of atmospheric gust velocity have 
been studied from airplane fiight measurements (ref. 2). 
Measurements of pitching velocity of a B-25 airplane in 
flight through rough air were used to determine the power- 
spectral-density function of atmospheric turbulence for four 
weather conditions. Clementson determined the reduced 
autocorrelation function R(t)/R(0) of pitching velocity and 
obtained the output power-spectral-density function by 
taking its Foimier transform. The appropriate frequency- 
response function as determined from simplified theoretical 
calculations was then used in equation (7) to obtain the 
input spectrum. The normalization of the autocorrelation 
function R{r) by dividing through by f?(0) is arbitrary and 
was made on the basis of yielding a pitching-velocity power 
of 1 (deg/sec) ^ for the airplane used in that investigation. 
The normalized power-spectral-density fimctions of atmos- 
pheric turbulence derived in this manner did not vary 
appreciably between weather conditions; although, as might 
bo expected, the value of i2(0) for the pitching velocity out- 
put did vary appreciably, and reflected variations in the 
average power of turbulence with weather condition. On 
the basis of the residts obtained in reference 2, the conclusion 
was reached that “atmospheric turbulence is a stationary 
random process that can be statistically described by a 
single reduced power-spectral-density cmwe.” Although 
additional tests under a wide variety of atmospheric condi- 
tions and the use of other airplanes are needed to verify this 
conclusion, the spectrum obtained appears to be representa- 
tive at least of the conditions covered by the tests. The 
results obtained in reference 2 for the power-spectral-density 
ftmction of atmospheric turbulence thus provide a tmbulence 
input and are used in subsequent applications in the present 
study. 

The average reduced power-spectral-density function ob- 
tained in reference 2 for an airplane true flight speed of 300 
feet per second is shown in figure 1. The results shown are 
corrected for two errors made in reference 2 and pointed out 

308CC5— 00 20 
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Angulor frequency, ii>, rodions/sec 


Fioxtbe 1. — Normalized power-spectral-denBity function for atmos- 
pheric vertical gust velocity 4>j(£i)) for airspeed of 300 feet per second. 


in reference 7. In addition, the results of reference 2 have 
been divided by 2ir ia order to conform to the definition of 
the power-spectral-density function used herein, equation (6). 
The power-spectral-density fimction shown in figure 1 has 
(ft/sec)^ 


dimensions of 


radians/sec 
speed of 300 feet per second. 
Changes in the frequency 


and applies to an airplane air- 


argument . — As pointed out 


previously, $(o)) has the dimension of a; 


power 


■and thus 


' radian/sec 

depends upon the airplane forward speed. In order to ex- 
press the power-spectral-density functions in terms inde- 
pendent of airplane flight speed, the change of variables 


«=7 


x=Vt) 


( 11 ) 


is introduced. The variable Q is a reduced frequency in 
radians per foot. The variable x is the airplane flight dis- 
tance in feet. In terms of these variables, the average 
power of a disturbance ^(x) is given by 

lim 4 r {^{x)fdx (12) 

X-*<o c/0 

(The use of the circumflex is explained in the section 
“Symbols.”) 

If the power-spectral-density function $(J2) is given by 


$(n)=lim f ^(x)e~*°^dx 

X-*co TT-A. ^0 


(13) 


then $(fl) gives the average power of the disturbance arising 
from components having frequencies between Q and Q-\-d^l. 
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J $(Q)cto=J ^(co)do}=^ 


it follows from equation (11) that 

l(Q)=y$(w) = F$(FQ) 


(14) 


(15) 


As pointed out in reference 7, the power-spectral-density 
fimction for atmospheric tm-hidence representing $j(Q) is 
incorrectly given in reference 2. The power-spectral-density 
function $i(Q), which corresponds to the function shown iu' 
figure 1, was obtained by equation (15) and is shown ia figure 
2. The resulting power spectrum of atmospheric turbulence 
is iadependent of airplane speed and in a limited sense may 
represent a universal reduced power-spectral-density function 
of atmospheric turbidence. 

The relation between the frequency-response functions for 
the sinusoidal input e*“‘ and is now considered. Since for 
a given airspeed the gusts represented by e*“‘ and are, 
from equations (11), the same gusts, the load responses when 
expressed in tencs of time t and distance x wfil only involve 

a change of scale. If !?'(if2)e'“* is the airplane load response 
to the gust then the frequency-response functions are 
simply related in the following manner: 

^'(iQ)=T('U>)=T(iVQ) (16) 

The frequency-response fvmction may be obtained from 
relations similar to those of equation (S) and given hy 

i*(tQ)=J^”-^(a:)e-‘“*da: 

!f^(iQ)=faJ' A(x)e~*^dx 

where tt^(a;) and A(,x) are the response functions of the unit 
impulse and unit step disturbances, respectively, expressed 
as a function of fiight distance. These response functions 
are related to W(t) and A{t), respectively, hy the relations 

^(x)=^ W(f) 



A(x)=A{t) 

In terms of the variable 0, the input-output relation can, 
from equations (15) and (16), be expressed by the relation 

$„(0)=$,(Ji)17’(iQ)P (18) 

The average power for a given input and output in terms of 
the frequency Q may from equations (14) and (18) be ob- 
tained by the relations 

^2(x)=J^ $<,(n)dQ=J^ $^(Q)i;?‘(ifi)Pdil (19) 

Another set of variables is useful in the applications of 
power-spectral analysis to gust loads. This set involved the 
reduced frequency v in radians per chord and the nondimen- 
sional distance s in chords. Because of the need for con- 
sideration of unsteady acrodynarrics in gust-loads analysis. 



Fioubb 2. — Nonnalined powei^spectral-density funotion for atmos- 
pheric vertical gust velocity 


calculated responses to step gust velocity inputs may fre- 
quently be given in terms of the nondimensional distance a 
and the frequency-response fimction in terms of the reduced 
frequency v. These variables v and a ore related to w and il 
and to t and x in the following manner: 



uc ^ 

y— 1^ = 011 


( 20 ) 


'If T{iv)e*'* is the acceleration response to the gust e'", then 

T'(iv) may also be obtained from the impulse or stop response 
by the relations 


T{iv)—^ W(a)e *'‘da 

T('ij')=ir f A(a)e~*'^da^ 
Jo , 


( 21 ) 


where lT(s) and A(s) are the responses to the impulse and 
step gust inputs in terms o'f the variable a. These response 
functions are related to W(f) and A{t), respectively, by the 
relations 


W(s)=^ TF(0 

A{a)=A{t) 


The frequency-response function T{iv) is related to the 
functions T{jA) and !f*(fO) in the following manner: 

^(tfi) = T(icG) = T(iFfi) (22) 

and involves only a change of scale of the frequency axis. 
These relations are used in the subsequent applications. 


EELATION OF POWEE-SPECTHAL-DBNSITY FUNCTION TO APPLIED GUST 

LOADS 

In order to apply the foregoing methods of analysis to the 
study of gust loads, the power spectral density of loads must 
be related to the intensities of the actual loads. The power- 
spectral-density function, which provides a measure of the 
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average power arising from components at various frequen- 
cies, does not directly reflect the load intensities since the 
actual load at a given time represents the combined output 
at the various frequencies. Thus, it is desirable to relate 
the power spectrum of loads to specific quantities of con- 
cern in load studies, such as the proportion of total time at 
a given load intensity (probability distribution of load in- 
tensity), the number and intensity of peak loads, and other 
such particular quantities that may be of interest for struc- 
tural analysis. In the present report, only the relation of 
the power-spectral-density fimction to the probability 
distribution of load is considered in detail and more specifi- 
cally the significance of the normal distribution for loads 
is considered. 

Probability distribution of output. — When a linear system is 
exposed to an input varying m a random manner with time, 
the probability distribution of the system output y can 
frequently be represented by a normal probability density 
distribution defined by the relation 



(23) 


the power spectrum alone defines a number of useful statis- 
tical characteristics of the time history.’ Some of these 
Telations between the spectrum and the time history statis- 
tics are discussed in the following paragraphs. 

Significance of the normality of the distribution of loads. — 
If the distribution of airplane acceleration mcrement in 
continuous rough air is normal with a zero mean value, the 
probability density distribution of the acceleration incre- 
ment An is completely described by the standard deviation 
and from equation (23) is given by 





(26) 


Equation (25) is termed the probability density fimction of 
An, and /(An) d(An) can be considered to represent the pro- 
portion of total time (or total flight distance) that An has 
a value between An and An-l-d(An). 

In order to examine the properties of equation (25), it is 
convenient to consider the distribution of the variable 



( 26 ) 


where y and a are the mean and standard deviation and are 
defined by 


where z is the so-called standardized variable and has the 
probability density distribution 


y= Urn i f y dt (24a) 

T-.® 1 Jo 

Investigations of communication problems associated with 
noise, which has many obvious similarities to turbulence, 
have shown that normal distributions are frequently en- 
coimtered. Rice in reference 8 has, for example, shown that 
for a linear system the shot effect in a vacuum tube gives 
rise to a noise current which has a normal distribution of 
current intensity. Investigation of fluid turbulence fre- 
quently yields normal distributions of velocity fluctuations. 
These results appear to be explained by the central limit 
theorem of probability (ref. 9) which states that, under gen- 
eral conditions, the distribution of the sum of a large num- 
ber of random variables tends toward a normal distribution. 

The theoretical derivation of the probability distribution 
of loads for an airplane in flight through rough air is math- 
ematically difficult and appears to involve assumptions 
regarding the nature of turbulence that are, for the present, 
questionable. However, as previously noted, turbulent 
velocity fluctuations have been observed to approximate 
a normal distribution. For the condition of normally dis- 
tributed velocity fluctuations, it may be expected that the 
loads experienced by an airplane in flight through the 
turbulent air would likewise be normally distributed. Some 
expeiimental data on the normality of the frequency dis- 
tribution of loads will be presented in a subsequent section 
of this report. The significance of the question of normality 
stems from the fact that, for the case of a normal distribution. 


^(s)=-^e-**/^ (27) 

V2ir 

The function ^{z) is the normal distribution with a mean of 
zero and a standard deviation equal to 1 and is the error 
function commonly tabulated. The probability that a 
random value of z will exceed a given value is given by the 
integral of equation (27) as follows: 

P(z)=J^ ^{z) dz (28) 


Equation (28) defines the cumulative probabflity distribu- 
tion of the reduced variable z and in the case of time-history 
data may be considered to represent the proportion of total 
time that the value of z exceeds a given value. For a fixed 
value of probability P, z is fixed and can be obtained from 
tabulations of the integral of the error function. Thus, for 
example, for P=0.02276, z=2. 

Equation (28) indicates that the probability of exceeding 
a given value of 2 is a function of only the given value. Con- 
versely for a given value of probability, the largest value of 
z exceeded is also fixed and depends only upon the probabil- 

ity. Making the substitution z = — into equation (28) 

yields the result that the probability of exceeding a given 
value of An/(T is likewise a function only of the value of Aw/o- 


and is given 



The converse also applies that, for 


a given value of probability, the largest value of An/o- ex- 
ceeded is fixed and depends upon the probability. Con- 
sequently, the largest value of An exceeded with a given 
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probability depends upon the value of o- and is given by <rz. 
The largest value of An exceeded with a given probability, 
is thus seen to be directly proportional to the standard de- 
viation. For example, when the standard deviation is 
doubled, the largest value of An exceeded with a given prob- 
ability is also doubled. This direct relationship between . 
the value of the standard deviation and the load probability 
distribution makes the standard deviation an important 
and significant measure of the load experience for the case 
of gust loads' having a normal distribution. Because of this 
direct relation between the standard deviation and the load 
probabUity distribution, the standard deviation will be useful 
as a measure of load intensity in the present study. 

The foregoing discussion has served to establish the plausi- 
bility of gust loads having a normal distribution under some 
conditions and the significance of the standard deviation of 
loads for the case of a normal distribution. The standard 
deviation of loads can be derived from a power-spectral 
analysis in the following manner: From its definition (eq. 
(24b)), the standard deviation of load increment output 
having a zero mean value is the square root of the average 
power. Thus, the power-spectral-density fimction of loads 
and the standard deviation of the probability distribution of 
loads are related. The standard deviation cr may thus be 
obt-ained directly from the power-spectral-density functions 
by the relation 

(T*= r $o(w)dw= f 4(Q)dQ (29) 

Jo Jo 

This relation between the probability distribution and power- 
spectral-density fimction for the case of a normally distrib- 
uted output ties the power spectrum to a basic characteristic 
of the load history and is thus of importance for applications 
to gust-load analysis. 

In addition to the foregoing relation between the power 
spectrum and the probabfiity distribution, a number of other 
useful relations exist for particular cases between the power 
spectrum and the statistical characteristics of the time 
history. Rice in reference 8 has derived expressions which 
permit the determination from the power spectrum of such 
quantities as the average number of crossings per second of . 
given values of disturbance intensity and the average 
number of peak values per second exceeding given values of 
intensity. These quantities are of particular interest in 
studies of airplane fatigue fife. These additional relations i 
substantially increase the usefuUness of the power spectrum. 

APPLICATION TO GUST-LOAD PROBLEMS 

In view of the simple relations between the loads and the 
power-spectral-density functions for the case of normal 
distribution of loads, the determination of the normahty of 
load distributions appears to be an important problem in the 
appUcation of power-spectral analysis to gust-load problems; 
therefore, some experimental gust-load distributions are 
examined for normality. Two applications of power-spectral 
methods to the calculation of gust loads are also given in 
this section in order to obtain some insight into the relations 
between airplane characteristics and loads in continuous 
rough air. In the first series of calculations, the input- 


output relations are used to calculate the standard deviation 
of loads for a selected series of idealized responses to a step 
gnst. The second series of calculations is made to determine 
the variations in the standard deviation of load for indi- 
vidual variations of certain geometric and aerodynamic 
parameters of an idealized transport airplane. The indicated 
variations in load intensity in this application are compared 
with those derived by conventional techniques of using the 
peak load response to an idealized discrete gust. 

OBSERVE® DISTRIBUTIONS OF LOADS 

In order to determine whether frequency distributions of 
'load are actually normal distributions, recourse was made to 
some available experimental gust-load time-history data. 
Time-history records of the normal acceleration for two au- 
planes of the same type (differing only shghtly in center-of- 
gravity position) in side-by-side flight through continuous 
rough air at low altitudes above generally flat terrain were 
available from a recent investigation. The test conditions 
and summaries of the data are presented in reference 10. 
The normal-acceleration time-history, records for about a 
2-miaute section of one run at 460 miles per hour were evalu- 
ated in detail by taking readings at intervals of second, 
roughly one reading for each 4 chords of airplane travel. 
The frequency distributions for the two airplanes are sum- 
marized in table I. (Additional data of the same type were 
also examined but do not appear to warrant reporting in 
detail at this time.) In determining these distributions, the 
1 g level was assumed to be at the mean value of the disturb- 
ance. (This assumption is frequently used in gust evalua- 
tions in view of the diflBculty in exactly determining the 1 g 
line for the flight condition.) The other primary charac- 
teristics of the observed frequency distribution, standard 
deviation < 7 , coefl&cient of skewness « 3 , and coefficient of 


TABLE I.— FREQUENCY DISTRIBUTIONS OF 
ACCELERATION INCREMENT 


Airplane A 

Airplane B 

Acceleration increment, 
An, g units 

Number 

Acceleration Increment, 
An, g nnits 

Number 

—a 777 to —a 727 

2 

-aTlOto -aC69 

1 

—.727 to —.677 

0 

-.660 to 

-. 019 

1 

-.677 to -.627 

1 

-.619 to 

-.669 

2 

-.627 to -.677 

0 

-.669 to 

-.619 

4 

-.677 to -.627 

2 

-. 619 to 

-.469 

11 

-.627 to -.477 

2 

-.469 to 

-.419 

12 

-.477 to -.427 

7 

-.419 to 

-.309 

22 

-.427 to -.377 

11 

-.369 to 

-.319 

36 

-.377to -.327 

10 

-. 319 to 

-.269 

00 

-.327to -.277 

33 

-.269 to 

-.219 

99 

-.277to -.227 

67 

-. 219 to 

-. 169 

no 

-.227tO -.177 

101 

-. 169 to 

-.119 

200 

-.177 to -.127 

167 

-.119 to 

-.069 

216 

-.127 to -.077 

211 

-.069 to 

-.019 

264 

-.onto -.027 

316 

-. 019 to 

.031 

260 

-.027to .023 

386 

.031 to 

.081 

277 

.023to .073 

338 

.081 to 

.131 

216 

.073to .123 

240 

.131 to 

.181 

172 

.123tO .173 

140 

.181 to 

.231 

126 

.173to .223 

125 

.231 to 

.281 

77 

.223to .273 

65 

.281 to 

.331 

60 

' .273to .323 

43 

.331 to 

.381 

30 

.323to .373 

16 

.381 to 

.431 

22 

.373to .423 

10 

.431 to 

.481 

16 

.423to .473 

0 

.481 to 

.631 

0 

.473tO .623 

7 

.631 to 

.681 

2 

.623to .673 

1 

.681 to 

.031 

2 

,673to .623 

1 




Total 

2339 

Total 

2341 

An 

0 

nn 


0 

O’ 

01606 

9 


0.1803 

ca 

-0049 

oa 


-0. 013 

a\ 

4.^ 

a| 


3.360 
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kurtosis ai are also given in table I where the statistical 
characteristics of the distribution were determined by the 
follo'wing relations: 


(30a) 

1 s(Aa— Aa)® 
N 

(30b) 

1 s(Aa-A»)^ 

(30c) 


where N is the number of observations. For a normal dis- 
tribution tt 3=0 and «4=3. The distributions of table I are 
shown as frequency polygons in figures 3(a) and 3 (b). Also 
shown in the figures are the normal distributions fitted in 
accordance ■with standard statistical procedures by means 
of the calculated standard deviations of table I. 

The fit of the normal distribution to the data for airplane A 
sho'ws some tendency of the normal distribution to under- 
estimate the concentration of values about the mean of the 
distribution, whereas the fit of the data of airplane B appears 
generally good. In order to test the hypothesis that the 
observed samples are from a normal distribution, statistical 
procedures were applied to the observed statisticalparameters. 
The results of these tests indicated that the distribution of 
An for airplane B could be considered from a normal dis- 
tribution, whereas the distribution for airplane A would not 
be assumed from a normal distribution. The magnitude of 
the departure from normality did not, however, appear large 
in this case. 

In order to permit examination of the beha'vior of the 
distributions at the larger values of load factor which are of 
concern, the observed relative cumulative frequency dis- 
tributions and the fitted cumulative probability distributions 
are shown in figure 4. The curves for both the probability 
of exceeding given values and the probability of being less 
than given values are sho-wn on semilogarithmic paper in 
order to permit comparison at both the large positive and 
negative acceleration increments. Examination of these 
figures indicates that the data for airplane B are in excellent 
agreement with the fitted curves. For airplane A, the 
overall agreement bet-ween observed data and fitted cmrve 
appears reasonable although some discrepancy between the 
observed distribution and the fitted curve is apparent par- 
ticularly at the laiger negative values of An. These dis- 
crepancies are, in general, not large, however, and might, 
for example, be due to piloting-technique effects. However, 
the lack of consistency in the results indicates the need for 
further study of the question of the normality of the dis- 
tribution of load. 

In order to examine the question of normality somewhat 
f'urther at this time, recourse was made to additional data 
on frequency distributions of load increment 'under condi- 
tions similar to those represented in figures 3 and 4. Exami- 
nation of these additional data indicates that although the 
distributions were generally close to normal, in a number of 
cases, more large accelerations were experienced than might 
be expected for a normal distribution. Consideration of the 


time-history records from •which these distributions ■were 
obtained indicated that these departures from normality 
were associated with lack of homogeneity in the turbulence 
intensity dming the flight run. Only for short runs of 
constant level turbulence did the distributions appear normal. 

The manner in which the departures from normality arise 
in nonhomogeneous turbulence may be clarified by an ex- 
ample. Consider a flight run through rough air consisting 
of somewhat more severe tobulent conditions during the 
second part of the nm. The overall distribution of loads 
may for this condition be expected to consist of two normal 
distributions, one for each part of the run. The two normal 
distributions would be expected to have mean values of zero 



(a) Airplane A. 

(b) Airplane B. 

Figobb 3. — Comparison of observed frequency distribution with fitted 
normal frequency distribution. 
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for tihe load factor increment but different standard devia- 
tions. The combined distribution for the whole run can be 
shown (appendix) to depart from normality with an excess 
of observations at the center and at the laager values of 
load increment. This mechanism can accoimt for the de- 
partures from normality observed in. the data examined. 

In view of the foregoing indications, it would appear 
reasonable for the present to assume that for locally homo- 
geneous turbulence of the type represented in figures 3 and 
4 the frequency distribution of load increments may be 



(a) Airplane A. 

(b) Airplane B. 

Fiqtjeb 4. — Ck)mpariBon of observed relative cumulative frequency 
distribution vritb fitted normal probability distribution. 


expected to approximate a normal distribution. The approx- 
imate normality of the load distribution for homogeneous 
turbulence permits the use of the relations derived earlier 
between the frequency distribution of loads and the power 
spectrums and permits the representation of the load inten- 
sity by a single munber, the standard deviation of the fre- 
quency distribution of load. This unification of the powor- 
spectral-density function and the probability distribution of 
loads even for limited conditions is of considerable importance 
since in many load studies the relative loads in continuous 
roigh air of two airplanes (or one airplane at two flight 
conditions) are of interest. For this case the use of the 
standard deviation as derived from the power-spectral analy- 
sis appears to provide a direct measure of the relative loads. 

In view of the foregoing indications that a lack of turbu- 
lence homogeneity can cause significant departures from a 
normal distribution of loads, further studies of the distribu- 
tion of loads under various atmospheric conditions are needed. 
It may be expected that the turbulence connected “with such 
dynamic phenomena as thunderstorms which have a large- 
scale physical structure and short life cycle may not be 
stationary random. For these conditions, the airplane load 
histories might have distributions that depart significantly 
from normality. This limitation may not, however, be too 
serious since even in these cases the standard deviation of 
loads may stiU provide a measure of load intensity although 
one not as simply interpretable in terms of the probability 
distribution as in the case of the normal distribution. 

EKLATION OF STEP EESPONSB CHAEACTEEISTICS TO LbADS IN 
CONTINUOUS BOUGH AIB 

Examination of calculated responses to entry of a sharp- 
edge gust for airplanes representative of conventional trans- 
port types has indicated that the shape of the response 
curve up to the peak acceleration depends primarily on the 
gust penetration function (Kilssner function). After the 
peak value which occurs dose to 6 chords of gust penetra- 
tion, the character of the response appears to be primarily 
a function of airplane stabihty and to approximate the short- 
period oscillation of the airplane. On the basis of these 
properties, a limited series of response curves to a sharp-edge 
gust were sdected to represent variations in airplane dy- 
namic longitudinal stability. The incremental acceleration 
responses in p’s to a 1-foot-per-second sharp-edge gust were 
assumed to be given by the following expressions: 

An(s)=isin^s 

An{8)=-^ e~’’ <*“®’cos io(s— 6) 

Equations (31) represent a quarter sine wave up to a fixed 
peak value of 0.033^ at 6 chords and a damped oscillatory 
function for the remaining portion of the response. The 
parameters 6 and can be considered to represent the 
damping and frequency parameters, respectively, of the 
response following the peak load. In figure 5, plots of equa- 
tions (31) are shown for nine selected cases covering three 
values for each of the two parameters. The values were 
sdected to sample a wide range of airplane response charac- 
teristics and represent variations in wave length from 40 to 
150 chords and variations in damping from light to almost 
critical damping. 



(31) 


Ln, g units 
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(a) *0=0.04 

(b) *0=0.07. 

(c) *0=0.10. 

FiauEB 5. — Load-faotoivincrement responses due to penetration of a sharp-edge gust for several stability conditions. 
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In view of the linearity of the systems hemg considered, 
the multiplication of the response An{s) of equations (31) by 
a constant results in the multiphcation of the frequency- 
response function aud the standard deviation of the output 
in turn by the same constant. Thus, the effect on the 
standard deviation of load of variations in the peak values 
of An for given values of ko and b is apparent and the con- 
sideration of an amplitude factor in equations (31) does not 
warrant more detailed discussion. 

The frequency response function T(iv) was determined 
analytically from the sharp-edge-gust response given by 
equations (31) by using equations (21). The amphtude 
squared of the transfer function ir(iv)l* is given by 





6vsin6v-f(6“— io^+r^r^cosev* . 


1 )(i*-l-io*+r^ 6 i<co86v— ( 6®— to*-t-v^j^sin6r 



Plots of equation (32) are shown in figure 6 for the nine 
cases considered in figure 5. 

The reduced power-spectral-density fimction of atmos- 
pheric turbulence, shown in figure 2, was used as the input. 
For any particular weather condition, the power-spectral- 
density fimction would have to he adjusted by t ak i n g into 
account the actual average power of the input. This change 
would, however, only involve an appropriate linear change in 
the ordinate scale and would not affect the relative indications 
of the present comparisons. 

Examination of figure 6 indicates that, as the damping of 
the oscillation is decreased (b decreased), the amplitude 
increases rapidly in the neighborhood of the oscUlation 
frequency ko. The frequency response function remains 
relatively unchanged over the rest of the frequency range. 
On the other hand, variations in the oscillation frequency ko, 
for fixed values of b, have a minor effect on the shape of the 
function but the peak values change (note changes in the 
ordinate scale) and occur at values of v dose to ko- 

The power-spectral-density function of the accderation- 
iucrement output was obtained from equations (18) and (22) 
by the relation 

$,(Q)=^t(fi)|T(wO)|* (33) 

For convenience, an airplane chord of 9.67 feet, the same 
value as the chord of the airplane used to obtain the iuput 
spectrum, was assumed for this series of calculations.. The 


power-spectral densities of the acceleration-increment output 

$«(Q) for the conditions being considered are shown in figure 7. 
Since the power-spectral-density function of atmospheric 
turbulence was not known for frequencies less than 0.0016 
radian per foot, the output spectrum could not be determined 
in this region. The output spectrums were, however, 
extrapolated to zero at fl=0 in order to complete the output 
spectrum at the low frequencies. The extrapolations are 
indicated in figure 7. 

Examination of figure 7 indicates that as damping is 
decreased the output spectmms for each value of kg increase 
rapidly in the neighborhood of the oscillatory frequency 
k k 

ft=-2=— For given values of damping, however, the 

power spectrums do not vary appreciably in shape but the 
peak values shift in frequency with ko- Thus, the total power 
of the spectrum as measured by the integral of the spectral 
function appears to be largely independent of ko and piimarily 
a. function of b for the conations investigated. This result 
can be seen from figures 6 and 7 to be a consequence of both 
the variations in response functions and the rapid decrease 
of power of the gust spectrum with increasing frequency. 

The rapid increase of the loads with decreased damping is 
more clearly illustrated in figure 8, which shows the standard 
deviation of loads <r (the square root of the integral of the 



(a) fro=0.04. 

(b) *0=0.07. 

(o) *0=0.10. 

Fiqueb 6. — Modulus squared of transfer functions for selected varia- 
tions of the frequency parameter *o and the damping parameter b. 
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power spectrum) as a function of frequency for the several 
values of b considered. The standard deviations of the 
accelemtion increment a- were obtained from the output 
spectrums by using equation (29). In the evaluation of 
equation (29), the output spectrums as extrapolated from a 
frequency of 0.0016 radian per foot to 0 (fig. 7) was used and 
the area of the output spectrums for values of Q greater than 
0.028 was assumed negligible. The figure indicates dearly 
that the load level, as measured by the standard deviation, 
is largely independent of the frequency parameter ko but 
varies appreciably when the damping parameter b is changed, 
the variation of the standard deviation being of the order of 
two to one for the conditions considered. 

The results imply that the short-period response character- 
istics and particularly the damping characteiistics have an 
appreciable effect on the airplane loads in continuous rough 
air. In contrast, the airplane peak-load response to discrete 
gusts can be shown for the present example to be affected only 
to a minor extent by the short-period damping characteristics. 
In view of these differences between the indications of power- 
spectral and discrete-gust calculations, discrete-gust calcula- 
tions may not adequately reflect differences in the gust loads 



(a) to =0.04. 

(b) to =0.07. 

(c) to=0.10. 

Figure 7. — Poweivspectral-density function of load-factor output. 
3080DG— 00 26 


in continuous rough air between airplanes differing in 
dynamic stability characteristics. In particular these results 
indicate that discrete-gust calculations may not he adequate 
for the determination of loads in continuous rough air for 
modem high-speed airplanes which in contrast to the air- 
planes of the past have relatively poor short-period damping. 
More complete accounting of airplane short-period charac- 
teristics in gust-load analysis would appear to he required in 
these cases. 

The variations of the response to a step gust assumed in 
the present illustration represent idealized conditions in 
which the characteristics of the airplane response to a step 
gust were changed in a simple manner. In practice, the 
change of almost any airplane parameter will modify the 
response to a step gust in a complex manner. Thus, a 
change of an airplane parameter will affect the peak-load 
value, the location of peak, and both the frequency and 
damping of the subsequent oscillation. Consequently, the 
problem of optimum design for gust-load reduction is 
extremely complicated and beyond the scope of the present 
study. The effects of variations in airplane geometry on 
gust loads are, to some extent, indicated in the second 
illustration in which the complete changes in the response 
to a step gust for limited variations in each of selected 
parameters of an idealized airplane are considered. 

EFFECT ON LOADS OF SOME VAEIATIONS IN AIRPLANE GEOMETRIC AND 
AERODYNAMIC PARAMETERS 

In order to obtain an indication of the relation between 
some airplane geometric and aerodynamic parameters and 
loads in continuous rough air, the power-spectral-density 
functions of load based on the turbulence spectrum of 
figure 2 were calculated for selected variations in airplane 
characteristics. An idealized transport airplane was used 
as a basic configuration. Characteristics of this basic con- 
figuration are listed in table II. In addition, values of each 
of 11 airplane parameters listed were varied separately in 
order to obtain an “increased” condition and a “decreased” 
condition, as indicated in the table. (In the case of the 
center-of-gravity position, the increased condition represents 



Frequency porometer, Aq 


Figure 8. — Standard deviation of loads as a function of the frequency 
parameter kt for various values of the damping parameter b. 
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TABLE n.— VALUES OF AIRPLANE PARAMETERS 


[U=l fl/sec; F-303 tt/seo: ci-7^ ft] 


Airplane parameter 

Basic 

Increased 

(•) 

Decreased 

(•) 


6.56 

6.40 

4.73 


3.21 

3.60 

2.73 


738 

886 

600 


276 

330 

220 


a002049 

a002378 

a 001766 


0.5 

ae 

0.4 

Z ' ff 

0 

-L26 

**L28 

Z./ft 

36.38 

39.26 

31.54 

T, RhipJtl 

209,600 

10 

272^000 

1L6 

146,100 

8.6 


iV; ih 

38,000 

47,600 

28,600 



* Values Indicate parameter changed for given condition, whereas other parameters remain 
samo os basic condition, 
b Center of gravity rearward. 


the forward center-of-gravity position.) The 23 conditions 
covered in these calculations represented stable airplanes 
with static margins varying from —0.25c to —0.50c where c 
is mean aerodynamic chord. 

The responses to a unit step gust with a velocity of 1 foot 
per second .4(a;) for the two-degree-of-freedom case, vertical 
motion and pitch, were available for the 23 cases of table II 
from an impublished study based on the methods of 
reference 6. The frequency-response functions were deter- 
mined by means of equations (17). The power-spectral- 
density function for acceleration output corresponding to 
the turbulence input of figure 2 was obtained for each 
condition by use of equation (18). The standard deviation 
of acceleration increment was determined for each case 
from equation (29) by the relation 

v=[j;i(u)dQj" 

In performing this integration, it was assumed that the 
output for frequencies of 0 greater than 0.028 was negligible. 
The output spectrums were also faired to a value of zero 
at Q equal to zero as in the previous appheation. 

As a basis for comparing the results of the indications of 
the spectral analysis with the indications of conventional 
types of analysis based on responses to single representative 
gusts, a triangular gust with a gradient- distance (distance 
from zero to peak value) of 10 chords was selected as a 
representative gust condition. This gust condition repre- 
sents an average gradient distance for the more severe gust 
loads and is frequently used for analysis purposes as a 
measure of the airplane loads in rough air. The peak load- 
factor value ATiaax for a 1-foot-per-second triangular gust 
having a gradient distance of 10 chords was calculated for 
the two-degree-of-freedom case for each of the 23 conditions 
considered. 

The results of the power-spectral and sin^e-gust calcula^ 
tions are compared in figure 9 for each of the 23 conditions. 
The abscissa in the figure represents the single-gust peak 
response and the ordinate, the power-spectral-density 
measure of load intensity. For ease in comparing the rela- 
tive changes, both single-gust and speotral-calculation results 
are shown in terms of the basic airplane response. Thus, the 


abscissa is given as An„axKAn„ax) i»ate and the ordinate, as 
<rl<rnuic- The line of equality — =— ^ is shoAvn as 

K^'^maxjlnuie ^ha$io 

a reference. If both methods of analysis gave the some per- 
centage change in load for a given change from the basic /lir- 
plane, the plotted point would fall on the line of equality. 

Examination of the results shown in figure 9 indicates 
that to a first approximation both the standard deviation of 
loads and the maximum load for a 10-chord-gradient triangu- 
lar gust show the same trends for variations of the airplane 
parameters considered. The largest changes in load level 
are associated with changes in airplane weight, wing area, air 
density, and slope of the wing lift emwe as might be expected 
from the sharp-edge-gust relation. Inspection of the figmo 
indicates that complex second-order dififerences that may be 
important exist in the indications given by the two measures 
of loads. As an example, for the variations of loads for 
changes of wing area the figure indicates that, for the 20- 
percent increase in wing area (represented by the square), 
the value of An^l{An^a^i^ic increases about 18 percent. 
The relative standard deviation txl<Tbanc, however, increases 
only about 7 percent. This difference appears to be a con- 
sequence of the greater sensitivity of the power-spectral 
analysis to the changes in airplane stability introduced by 
the increased wing area; the increase in wing area results in a 
small decrease in the frequency of the response to a step 
gust but , a more pronounced increase in the damping. The 
present results also suggest similar dififerences between the 
two analyses for changes in airplane weight and slope of the 
wing lift cmwe. 

The variations considered for the remaining parameters 
such as tail length, tail area, and tail slope of the lift curve, 
in general, yield min or variations in loads. The variations in 
An^liATi^ai) KutpioT these cases are however less than ± 2 per- 
cent, whereas the variations for the standard deviations are 
generally somewhat laiger, ±4 percent. In figure 9 the 12- 
percent-chord rearward movement of the center-of-gravity po- 
sition la yields about a 3^ercent increase in An„axK^n„ax)tmtie 
but a 9-percent increase in the relative standard deviation. 
In this case, the larger increase in the standard deviation 
appears largely associated with the movement of the peak 
of the frequency-response function to lower frequency and 
thus larger gust spectral, power. 

For the stable airplane configurations considered, the 
power-spectral analysis and the sin^e-gust peak-load re- 
sponse yield trends for changes in parameter values that are 
generally consistent. Dififerences that may be important 
are, however, discernible and appear associated with the 
greater dependence of the standard deviation of loads on the 
dynamic stability characteristics of the airplane. It would 
be expected from the first application that, for larger varia- 
tions in the damping characteristics of the short-period oscil- 
lation, the dififerences between the two types of analysis 
would be appreciably larger than obtained in the present 
illustration. 




Stondord deviation ratio, 
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Fiqubb 9. — Comparison of standard deviation of load with peak load response to a 10-chord-gradlent triangular gust. 

Flagged symbols denote below-bMio. parameters. 
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SUMMARY OF RESULTS 

The analysis of the apphcation of power-spectral methods 
of analysis to gust-load problems has indicated the foUowiag 
results: 

1. The apphcation of power-spectral methods of analysis 
to load calcidations provides a measure of load intensity for 
continuous rough air in terms of the standard deviation 
(root mean square) of the load output. 

2. The probabhity distribution of load intensity in homo- 
geneous rough air appears to approximate a normal 
distribution. 

3. For the case of the normally distributed output, the 
standard deviation of load completely describes the prob- 
abihty distribution of loads specifying the proportion of total 
time that various load values are exceeded. 

4. The apphcation of power-spectral relations to a selected 
series of systematic variations in the frequency and damping 
characteristics of the airplane load response to a step-gust 
input indicates that the damping characteristics of the indi- 
cial response are the primary characteristics in determining 
the loads in continuoiis rough air. This result appears of 
signihcance in regard to high-speed airplanes and missiles 
where the short-period longitudinal damping may be poor. 

5. Calculations for a limited series of conventional and 
stable airplane configurations indicate that the loads in 
continuous rough air for variations in individual airplane 
geometric and aerodynamic parameters are to a first 
approximation adequately reflected by the peak-load response 
to the arbitrary 10-chord triangular gust commonly used. 


However, differences are discernible between the indications 
of the contmuous-gust and discrete-gust calculations and 
appear largely associated with differences in the effects of 
stabihty changes on the loads for the two gust conditions. 

SUGGESTIONS FOR FUTURE RESEARCH 

The foregoing results appear to indicate that power- 
spectral methods are weU suited to the calculation of loads 
in continuous rough air. A number of problems exist and 
require further investigation. These include 

(1) The determination of the power-spectral-density fimc- 
tions of atmospheric turbulence for a wide range of atmos- 
pheric conditions 

(2) The determination of the conditions under which 
normal distributions of the load time history apply 

(3) The investigation of the effects of the nonuniformity 
of gusts across the airplane span on the apphcation of these 
techniques 

(4) The experimental verification of theoretical determi- 
nations of airplane mdicial responses over the range required 
to determine usable frequency-response functions 

(5) The investigation of the relations between the power 
spectrums and the frequency distribution of load-increment 
peaks and other quantities of interest in structural design. 

Langley Aeronautical Laboratory, 

National Advisory Committee fob Aeronautics, 
Langley Field, Va., September S9, 195B. 


APPENDIX 

DISTRIBUTION OF GUST LOADS IN NONHOMOGENEOUS ROUGH AIR 


If the two normal frequency distributions of load incre- 
ment with zero mean are defined by 

-L 

iNTi/i(An)= — j=e 2 
(TiVStt 

N,MAn)=-^e-^ 
o'2V2'a- 

where 

fi(An),f 3 (An) respective probabhity density "functions 
NijNi respective number of observations for each 

distribution 

<Ti,t73 standard deviations of respective distributions 

the combined distribution is defined by 

Nf(An)=Nifi{An)+N2f2{An) (A2) 

where N = Ni + N 2 . This distribution is examined for 
normahty. 

Since the two frequency distributions (Aw) andiV 2 / 2 (A») 
are each normal with mean zero, the combined distribution 
Nf^An) is also symmetrical about zero. Consequently, for 
the moments of the frequency distribution defined by 

iia= I (Aw)"/(Aw) d{An) (A3) 

V — CD 



where m is a positive number and designates the moment 
order, the following relations apply 


/x„=0 (m odd) 

n„r^0 (m even) 


(A4) 


For a normal distribution, the coeflEicient of kurtosis 


defined by 

y-A yA 
“1— rr 

(A6) 

has a fixed value 

yi 



«4=3 

(A6) 


In order to examine the normality of the distribution 
defined by equation (A2), consider only the value of «<, the 
coefficient of kurtosis, for that frequency distribution. 

By definition, the coefficient of kurtosis a 4 for the distribu- 
tion of equation (A2) is given by 


j (Aw)Wi/i(Aw)d(Aw)-l-J (Aw)'‘N'2/a(Aw)d(A7i)J 


(Aw)“ [iVi/i (Aw)-|-A^*/a (Aw)] d(An) 

ivT+iVa 
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which when simplified yields 


Oi= 


3 

Nt+N, 


(N 1 <T j^+iVa <r 3*) 


A^i tr <ra^* 

. N^+N3 ) 


(AS) 


From equation (A6), a necessary condition for the distri- 
bution Nf{An) to he normal is that 04 be equal to 3. From 
equation (AS) this condition is obviously true only when 

N,a,*+N3<T3*=J^^^ (.N,a,^+N3<T3r (A9) 

Expanding the term on the right and simplifying yields 

a*+a3^=2aiW (AlO) 

Equation (AlO) is thus a necessary condition for the distribu- 
tion defined by (A2) to be normal. However, from the 
inequality 

tt*-f6*]>2a6 (ay^b) (All) 
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